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FOREWORD 

A  batter  understanding  of  sound  propagation  in  a  complex  medium 
like  the  ocean  or  atmosphere  continues  to  be  a  goal  of  vital  interest 
to  the  Navy.  In  many  practical  applications  it  is  sufficient  to 
approximate  the  ocean  as  a  stationary  medium  in  which  the  acoustic 
properties  vary  only  with  depth.  The  theory  then  is  very  simple. 

In  n  more  realistic  model,  however,  this  theory  is  imbued  with  a 
high  degree  of  complexity  due  to  general  three-dimensional  inhomo- 
geneities  as  well  as  motion  of  the  medium. 

This  report  addresses  Itself  to  the  exact  theory  of  ray  acous¬ 
tics  in  such  a  general  medium.  Differential  equations  for  the 
acoustic  rays  are  derived  in  two  different  forms,  which  are  suitable 
for  implementation  on  a  digital  computer. 

This  work  was  done  in  the  Science  Research  Grcup  of  the  Warfare 
Analysis  Department  under  Independent  Research  Project  Number 
R360FR103/2101/RQLI01001.  The  date  of  completion  was  l  June  1970. 

ft  p*  «  i  ■  ■»  T***  • 
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Abstract 

Two  different  formulations  are  presented  in  terms  of 
ordinary  second-order  vector  differential  equations  for 
acoustic  rays  in  a  three-dimensional  medium  moving  with  an 
arbitrary  velocity  (either  BubBonic  or  supersonic) ,  and  having 
an  arbitrary  Index  of  refraction.  The  first  formulation  has 
the  arc  length  of  the  ray  as  independent  variable,  while  the 
aecond  one  is  given  in  terms  of  a  canonical  variable  and  is 
equivalent  to  Keller's  Hamiltonian  formulation  [J.  B.  Keller.. 

J.  Appl.  Phys.  2£,  938-947  (1954)].  An  explicit  Lagrangian  is 
ccfti.it  met  for  this  problem*  from  which  a  generalized  Fermat* s 

principle  it  derived.  An  error  in  an  earlier  publication 
[P.  Ugincius,  J.  Acoust.  Soc.  Amer.  37,  476-47°  (1965)]  is 
corrected,  and  it  is  shown  that  this  error  disappears  when 
the  velocity  of  the  medium  is  small  compared  to  the  speed  of 
sound.  In  that  case,  and  with  the  further  restriction  of  planar 
rays,  an  explicit  expression  is  derived  for  the  curvature  of  the 
rays  in  terms  of  the  velocity,  speed  of  sound,  and  their  gradients. 


1 


INTRODUCTION 


In  Ref.  1  one  of  the  authors  derived  «  vector  differential  equation 
for  the  acoustic  ray  paths  in  a  moving  inhomogeneous  medium  i.i  which 
the  velocity  field  and  index  of  refraction  could  be  arbitrary  £unc= 
cions  of  position.  Unfortunately  u a  error  was  made  in  chat  derivation 
eo  that  the  main  result  14,  !Uf,  1)  is  valid  only  for  those  cases 

in  which  the  velocity  of  tbu  medium  is  much  smaller  than  the  speed  of 
sound.  Ia  this  paper  wo  derive  the  correct  equation  for  arbitrary 
velocities,  both  supersonic  a*id  subsonic.  This  is  done  in  two  ways, 
in  part  I  by  the  methc-i  of  fte£,  1,  and  in  Tart  II  by  finding  the 
characteristics  of  tha  time*- independe.it  eikonal  equation.  The  second 
method  of  derivation  is  important  for  two  reasons:  1.  being  based 
on  the  rigorous  theory  of  partial  differential  equations  it  demon¬ 
strates  the  validity  of  the  "tricky”  method  used  in  ?«f.  1;  and, 
mere  important,  2.  it  leads  to  a  natural  (canonical)  Independent 
variable  t  (different  from  arc  length  or  time)  in  terms  of  which 
the  differential  equation  for  the  ray  paths  assumes  a  simpler  form. 

The  latter  method  ia  quite  analogous  and  equivalent  to  the  Hamilton¬ 
ian  formalism  used  by  Keller8,  Keller's  work,  however,  does  not 
finish  the  problem,  because  it  does  not  consider  the  second  set  of 
the  canonical  equations  for  the  generalized  momenta. 

In  part  III  we  construct  an  explicit  Lagrangian  for  this 
problem  and  derive  a  generalized  Fermat's  principle.  The  correct 
expression  of  Fermat's  principle  for  a  moving  inhomogeneous  medium 
had  evidently  not  been  known  until  now*. 
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In  pact  IV  the  general  equation*  are  specialised  to  planar  rays  In 
a  medium  with  small  velocity,  and  it  is  shown  that  the  results  of 
Ref,  1  ere  correct  in  this  case. 

I.  DERIVATION  GF  THE  DIFFERENTIAL  EQUATION  FOR  THE  RAYS 
Let  the  acoustical  properties  of  the  medium  be  described  by 
the  speed  of  sound  c(j$ ,  and  let  its  velocity  be  %(£) ,  where  both 
c  and  ^  are  arbitrary  wall-behaved  functions  of  the  position  vector 
£.  He  shall  use  the  following  terminology: 

i/  “  <fc/d»  ■  unit  tangent  to  a  ray. 
n  ■  c^/c  ■  Index  of  refraction. 

A*  2/c  *  nondimens  tonal  velocity, 
i  -  ■  unit  norms 1  to  the  wave  front. 

Here  s  is  the  arc  lrngth  measured  along  the  ray,  c0  is  an  arbitrary 
constant  reference  velocity,  end  ft  is  the  elkonal  function  (phase) 
of  the  wave  front. , 


is  well  known*"*,  and  can  be  derived  as  the  short-wavelength  limit 
(geometrical  approximation)  of  the  wave  equation: 

-  n-d  -  jx  •  jtft/co)* .  (i) 

This  can  be  written  in  terms  of  our  terminology  as 

jvsj  -  p  -  ±(n  -  ^  •  £)  .  (2) 
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We  define  now  the  direction  of  the  acoustic  ray  path  to  be 
the  vector  sum  of  the  unit  normal,  £  "  /  p  and  the  velocity^  (see 

Fig.  1): 

£ +  x  ■  »  O) 


where  q  will  be  called  the  (nondimenslonal)  ray  speed.  (Actually  the 
true  value  of  the  ray  speed  is  cq).  This  geometrical  definition  can 
be  shown1*  to  be  identical  with  the  physical  definition  that  the  ray 
path  be  the  direction  of  energy  transport.  Multiplying  Eq.  3  by  p, 
we  have 

-  *&'  -  X*  *  <4> 


where  we  have  Introduced  the  two  new  auxiliary  functions 

N  ■  pq  }  ■  pv,  .  (5) 


The  desired  differential  equation  for  Che  ray  paths  is  obtained  by 
taking  the  directional  derivative  of  Eq.  4  in  the  direction  of  , 
and  then  eliminating  {  between  these  two  equations0: 


b^  ’  V  *  •  *‘V‘Z5  •  • 


ds 


(6) 


Knf-  A  y.TA  hairs  4  nf  AffhanAa^  f-Vte  V7  sm/1  X7  ava  4  ^  Vn  £  J 

*•*•*•■  *ilwtsW»*iM»evw  biiw  *  J  ******  ^  Vt'VkPbVkC  All  O  ailUj 

therefore ,  is  subject  to  the  usual  continuity  restrictions. 

We  now  substitute  the  1-th  component  of  Eq.  4  into  the  last  term 
of  Eq.  b: 


d  dV 

t-(Nx')  -  -r*  -  x,'S(Nx/  -  V.)  -  VN  -  x/W. 
ut  —  as  ■  —  -  ■  * *-  • 


{7) 
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In  arriving  at  the  last  term  of  Eq.  7  we  made  use  oi  the  fats  that 
is  a  unit  vector:  x/x/  *  1  .  Equation  7  is  identical  with  an 
Intermediary  result  of  Ref.  1.  However,  at  this  point  an  error  was 
made  in  Ref.  1  by  taking  dj£/ds  -  r/»  ^  *a  *n 

explicit  function  of  both  the  position  and  the  ray  direction 
through  tha..  point,  so  that  the  correct  expression  for  the  directional 
derivative  of  V  must  be 

3  -  *  x-z'ix  .  (8> 

where  %  is  the  usual  gradient  operating  on  the  position  coordinates  xi 
only,  while  V.'  ■  jiid/dx/  operates  only  on  the  direction  coordinates 
x/  ,7  Equation  7  now  becomes 

S^>.'>  ■  +  -  s>  •  w 

The  last  two  teems  in  the  left-hand  side  of  Eq.  9  can  be  combined  into 
a  triple  vector  product,  which  yields  as  the  final  form  of  the  differ¬ 
ential  equation  for  the  rays 

7r0&'>  -  &*•?')£  +  x/  x  (V  X  £)  -  JN.  (10) 

Comparing  this  with  Eq.  14  of  Ref.  1  we  see  that  the  error  there  is  the 
neglect  of  the  term  (gf  •%')£•  Equation  10  is  an  ordinary  second-order 
vector  differential  equation  which  can  be  integrated  to  give  the  ray 
path  £( s)  as  a  function  of  its  arc  length  a  for  any  given  initial 
values  £(0)  and^^O), 
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D«  The  Functions  N  cr.d 


In  order  to  integrate  Eq.  10  we  must  express  the  functions  N  and 
^  (defined  in  Eq.  5)  explicitly  in  terms  of  £  and  Writing  Eq.  3 
as  ■  q j '  -  ^  and  squaring  yields  a  second-degree  algebraic  equation 
for  q  which  can  be  solved  to  give 

q  ■  ±  S  ;  (11) 

s  »  [l  -  va  +  <£'*jr)a]*  .  (12) 


Note  from  Fig.  1  that  the  geometrical  interpretation  of  S  is  the  pro¬ 
jection  of  the  unit  vector  on  the  ray  direction  S  ■  cos(a-!3)  . 

In  terms  of  the  absolute  velocities,  therefore,  cS  is  the  projection 
of  the  wave  velocity  cjJ  in  the  direction  of  the  ray.  The  requirement 
that  q  >  0  for  v  <  1  leads  from  the  geometry  of  Fig.  1  to  the  conclu¬ 
sion  that  only  the  upper  (+)  sign  in  Eq.  11  is  allowed.  Similarly 
for  v  *  1  one  can  show  that  the  upper  sign  must  be  used  whenever 

cos  a  ■  (£*^)/v  a  -  l/v  .  (13) 

If  the  Inequality  in  Eq,  13  is  reversed  then  the  lower  (-)  sign  should 
be  used  in  order  to  have  q  i  S,  It  is  shown  below,  however,  that  the 
reversal  of  this  inequality  is  strictly  forbidden  by  another  restric¬ 
tion  on  the  function  p,  so  that  the  upper  sign  in  Eq,  11  holds  under 
all  conditions?  This  additional  restriction  can  easily  be  shown  to  oe 
equivalent  to  the  physical  requirement  that  the  rays  be  confined 
within  the  Mach  cone 

cos  0  ■  *  (v3-l)*/v  .  (14) 


* 


b 


Equation  14  also  guarantee*  the  reality  of  the  function  S  defined  in 
Eq.  12. 

To  obtain  p^,  jj/ )  we  solve  Eq.  2: 

p  “  n/(^*j>  *  ^  »  (15) 

which  proves  the  restriction  in  Eq,  13,  because  ve  must  have  p  >  0. 
Substituting  for  ^  •  £  from  Eq.  3  (and  using  only  the  upper  signs 
from  now  on)  we  can  write  Eq.  15  as 

p  -  n/(l  -  va  +  %,»j)  -  n/[l  -  va  +  +  Qfc/*iJi)s3 

■  +  S)S  .  (16) 

This  can  be  written  in  a  more  convenient  form  after  multiplying 
both  numerator  and  denominator  by  -  S  and  using  Eq.  12: 

P  ■  n(S  -  v8)S  .  (17) 

The  function  N  of  Eq.  5  obtained  by  multiplication  of  Eqs.  11  and 
17  now  assumes  a  very  simple  form: 

N  -  n/S  .  (18) 
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C.  The  Final  Equation* 


The  functions  N  and  ^  are  now  expressed  entirely  in  terms  of 
£  anijj/  by  Eq's.  5,  11,  12,  17,  and  18,  so  that  the  differential 
equation  for  the  rays  (Eq.  10)  can  be  integrated  for  any  given 
velocity  and  index  of  refraction  n^j,).  For  easier  reference 

all  of  the  necessary  equations  are  collected  below: 


(19a) 

N  ■  n/S  ;  V  ■  ; 

(19b) 

S  -  [1  -  va  +  <£'&)*$  ; 

(19c) 

p  ■  n(S  -  ~  v3)S  . 

(19d) 

Note  that  for  v  >  1  and  for  a  ray  close  to  the  Mach  cone  (approach¬ 
ing  equality  in  Eq’s.  13  and  14),  S  will  approach  sero,  and  therefore 
both  N  and  will  grow  beyond  any  bound.  Also,  for  v  *  1,  independent 
of  how  close  the  ray  la  to  the  Mach  cone,  S  ■  ,  and  therefore  p 

and  ^become  indeterminate.  In  both  of  these  cases,  therefore,  one 
should  expect  difficulties  when  integrating  Eq.  19a  numerically. 

II.  THE  CHARACTERISTICS  OF  THE  EIK0KAL  EQUATION 
Courant  and  Hilbert®  have  shown  that  in  the  trivial  casejj£*0,  n«l 
the  sound  rays  are  the  characteristics  of  the  elkonal  equation,  Eq.  1. 
In  this  section  we  shall  prove  this  to  be  the  case  also  for  our 
general  problem.  Besides  showing  that  Eq.  19a  is  the  differential 
equation  for  the  characteristics  this  method  will  also  give  us  a 
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somewhat  simpler  equation  for  the  raya. 

The  theory  of  partial  differential  equations1  tells  us  that  any 


first-order  differential  equation  for  t 

H(X| ,  p, ,  «)  -  0  ;  pi  *  ,  (20) 

ie  solved  by  a  family  of  characteristics  which  are  curves  given  by 
the  following  system  of  ordinary  differential  equations: 

dx,/dT  •  SH/aPt  ,  (21a) 

dpj/dT  -  -  (aH/axj  +  pt3K/a»)  ,  (2ib) 

dS/dT  -  PjSH/apj  ,  (21c) 


where  T  is  an  Independent  parameter.  The  eilronal  equation  (Eq.  1) 
written  in  the  form  of  Eq.  20  becomes 

H  -  PjPj  -  (n  -  vjpj)*  -  0  ,  (22) 

and  its  characteristics  as  given  by  Eq1*,  21a  and  21b  are  easily  shown 
to  be: 

*t  ■  2[  Pj  +  (n  -  Pjv^v,  3  ,  (23a) 

Pj  ■  2(n  -  vJpJ)[(an/bx1)  -  PjOvj/axj)]  ,  (23b) 

where  from  now  on  a  dot  will  signify  differentiation  with  respect  to  t. 
Equation  21c  gives  the  value  of  the  eikonal  function  i  on  the  charac¬ 
teristics  which,  however,  is  not  needed  for  tracing  the  rays.  The 
significance  of  Eq.  21c  will  be  explored  further  in  the  Appendix. 
Equation  23c  is,  apart  from  a  scale  factor,  equivalent  to  Keller : 

Eq.  38.  Keller,  however,  does  not  consider  the  second  set  of  the 
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characteristic  equation*,  Eq.  23b,  which  era  needed  to  eliminate  the 
pt  from  Eq„  23a  in  order  to  ofcmin  a  differential  equation  for  the 
rays.  Equations  23  with  the  help  of  Eq.  22  can  be  written  more  con¬ 
veniently  in  vector  tonr.’ ~  • 


i-2^+pv)  , 

£‘  2p(^  -  Pj<V;)  . 


(24a) 

(24b) 


Comparing  5q.  24a  with  Eq.  3  we  6 se  indeed  that  the  direction  ^  of 
the  characteristics  coincides  with  the  direction  of  the  sound 
rays.  It  still  remains  to  be  shown,  however,  that  Eq'e.  24  are 
identical  to  Eq.  10.  To  do  thin  we  first  eliminate ^  between  Eq'a,  24 
by  solving  for  from  Eq.  24a: 


X  "  2**  “  %  *  (25) 

Taking  the  dot  product  of  this  equation  with  ^  and  using 

*  £  -  n  -  p  (26) 

from  Eq.  2,  allows  us  also  to  determine  explicitly  p  (the  magnitude 
of  the  vector  j>) : 

P  ■  (n  "  v*>  •  (27) 

4m 

Substituting  Eq.  25  into  Eq.  24b  we  get 

fcdi,-  P*)-  2K5"  -  Pj^O  -  (28) 

Equation  28  (together  with  Eq's,  25  and  27  which  are  needed  for  p  and 
P j  in  Eq.  28)  is  a  single  second-order  vector  differential  equation 
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for  the  ray  paths.  It  looks  somewhat  simpler  than  Eq,  19  and 
therefore  might  be  preferred  for  practical  calculations.  The  inte¬ 
gration  of  Eq.  28.  hcwevor,  will  give  the  rays  in  terns  of  the  as  yet 
undefined  canonical  parameter  t,  whereas  Eq.  19  gives  them  as  functions 
of  the  phyelcel  ere  length  s.  Whether  one  is  to  be  preferred  over  the 
other  Is  not  known  in  general,  and  might  depend  on  the  particular 
problem  under  consideration.  For  example,1'12  In  acme  situations  e 
constant  integration  stop  As  is  desirable,  in  which  case  the  formula¬ 
tion  of  Eqe.  19  would  be  preferred.  In  other  cases13,  however,  where 
the  ray  might  wind  Itself  into  a  point  of  singularity,  e  constant  As 
would  be  completely  inadequate  for  numerical  integrations. 

To  complete  the  identification  of  Eq.  28  with  Eq.  19,  and  at  the 
same  time  to  determine  the  physical  meaning  of  the  canonical  variable  t, 
we  make  e  transformation  of  the  independent  variable  from  t  to  i  in 
Eq.  28.  This  is  done  easiest  by  substituting 

&  -  £'(ds/dT)  (29) 

into  Eq.  25,  solving  for  ?  and  then  squaring.  The  result  can  be 
simplified  by  noting  that  jj^.'  |s  »  1,  and  by  using  Eqs.  26  and  27  to 
substitute  for  all  fictoxe  containing  p  or  •  A  little  algebra 
will  finally  yield 

ds/d.T  ■  2n/S  ,  (30) 

where  S  is  the  same  function  of  £  and  which  was  used  in  section  I 
and  is  defined  in  Eq.  19c.  Equations  29  and  30  can  now  be  used  to 
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rewrite  tsqc.  26  end  2/  with  s  ea  the  independent  variable: 

“(S'-nS'i)  .  <Jl*> 

Ww  >  w  »  »i  »  r 

p  -  n(S  -  -  v8)S  .  (31b) 

Note  that  Eq.  31b  egreea  exactly  with  Eq.  19d,  end  that  the  left  side 
of  Eq.  31a  agrees  with  the  left  aide  of  Eq.  7.  It  ia  straightforward 
but  rather  tedious  to  show  that  the  right  sides  of  Eqs.  31a  and  7 
agree  also.  Therefore  we  have  proven  our  contention  that  the  charac¬ 
teristics  of  the  elkonal  equation  are  tha  acoustic  rays,  and  as  a  by¬ 
product  have  obtained  an  alternative  (possibly  simpler)  formulation, 
Eqs.  26,  25  and  27,  for  tha  ray  differential  equation. 

It  is  interesting  to  note  that  part  of  the  troubles  anticipated 
with  a  numerical  integration  of  Bq.  19a  do  not  appear  in  the  new 
formulation  of  Eq.  28.  That  is,  for  v  >  1,  there  are  no  quantities 
which  blow  up  when  a  ray  approaches  the  Mach  cone.  For  v*l,  however, 
the  same  difficulty  persists  as  in  Eqs.  19,  because  p  (Eq.  27)  then 
becomes  infinite.  We  can  understand  better  now  why  the  first  diffi¬ 
culty  appears  In  Eqs.  19  and  is  removed  in  tha  formulation  of  Eq.  28: 
Eq.  30  shows  that  close  to  the  Mach  cone,  where  S  approaches  zero, 
ds/dT  approaches  Infinity,  so  that  an  inflniteseimal  integration 
step  dr  In  terms  of  tha  canonical  variable  corresponds  to  an  Infinite 

Swop  In  Cvnoi  who  gift  length • 
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As  was  mentioned  in  the  Introduction,  the  method  of  character- 

1-  m-  1  -X  J  .U  mm— mm  ..aA<l  *>  /-»  A  ITU  f"VlO  WAV  ■(  f  fl  1  &»•<  «  t"  ^  4  H 

BLILO  WlliUI  WOP  U8»»  ww  «*wa*vv  >  mv  •  **y  -  Cn«.w..v..  ... 

section  II  is  equivalent  to  the  Hamiltonian  method  used  by  Keller  . 
In  fact,  since  the  function  H  of  Eq.  22  does  not  depend  explicitly 
on  *,  the  differential  equations  for  the  characteristics,  Eqs.  21, 
are  precisely  Hamilton's  canonical  equations  of  motion  which  are 
well  known  in  classical  mechanics14.  The  function  HiXt.Pt), 
therefore,  can  be  regarded  as  the  Hamiltonian  of  some  mechanical 
system,  and  with  it  we  can  construct  the  Lagranglan  from  which  the 
ray  differential  equation  la  derivable  via  a  variational  principle 
(^grange's  equations). 

The  Lagrangi&n  L(xi,x()  is  given  in  terms  of  the  Hamiltonian 
H(*i.Pi>  by14 

L  »  XjPj  -  H 

■  *jPi  ~  P1P1  +  O  "  v,p,)a  .  (32) 


In  order  to  put  the  Lagrangian  in  its  proper  form  (a  function  of  the 

1  ip)  taa  SAliro  f*k  n  uwa  >  aw  a  f  Prts  91?  o  *•  fha  n.  <  rt 

J  »HW  A}  WV  MU«  w  wv*v  V  Wttv  WWWtA*  W  X-  t  Wa.w  |P| 

terms  of  the  xj  (without  using  the  "energy"  equation,  Eq,  22),  and 
substitute  the  result  into  Eq.  32.  The  algebra  is  rather  lengthy, 
so  we  shall  omit  it  here  and  just  give  the  final  result; 


/ae\ 


VJJ^ 
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Lagrange's  equations 


J  /vr  \ 


(34) 


can  be  used  to  derive  the  ray  differential  equation  (Eq.  28)  once 
more,  proving  that  Eq.  33  ie  indeed  the  correct  Lagrangian  for  this 
problem. 

Having  obtained  the  Lagrangian  we  can  now  examine  the  varia¬ 
tional  principle 

B 

6  J*  L  dT  -  0  ,  (35) 

A 

(whose  solution  is  Eq.  34)  which  for  the  case  of  a  motionless  medium 
gives  the  well-known  Fermat's  principle,  which  states  that  the  ray 
path  between  any  two  points  A,B  is  that  curve  along  which  the  time 
of  travel  is  minimum.  Korahauser®  implied  that  such  a  Fermat's 
principle  may  not  exist  for  a  moving  medium.  He  also  asserted  that 
Rayleigh's1®  contention,  that  the  velocity  to  be  used  in  Fermat's 
principle  is  the  speed  of  sound  plus  the  component  of  the  fluid 
velocity  in  the  ray  direction,  is  Incorrect.  To  see  whether  Eq.  35 
is  a  Fermat's  principle  (in  the  sense  of  minimizing  a  travel  time) 
we  rewrite  the  Lagrangier.  it  terms  of  the  arc  length  variable  s 
(using  Eqs.  29,  30  and  3 lb) 

L  ■  2np  ,  (36) 
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whereupon  Eq.  (35)  can  be  written 

B  B  v  » 

S  j‘  L  dr  -  6  J  2np^l)da  -  6  pSds 

B 

-  fl  f  - -  —  ds  -  0  ,  (37) 

where  in  the  left  atop  we  have  used  Eq.  (lb)  to  substitute  for  p. 

The  denominator  of  thla  laat  expression  is  just  the  ray  speed  q 

defined  in  Eq,  11  so  that  the  variational  integral 

B  J 

6  f  —  ■  0  (38) 

A  C1 

is  indeed  Fermat's  principle  in  the  strict  sense,  because  it  mini¬ 
mizes  the  ray's  travel  time.  Referring  to  Fig.  1  and  using  Eq.  3 
we  see  that  the  velocity  to  be  used  in  tne  denominator  of  Eq.  38  is 
cq  ■  c  cos(or  -  8)  +  v  cos  8,  whereas  Rayleigh' sXB  idea  was  to  use 
c  +  v  cos  8.  Kornhauser®,  therefore,  was  correct  in  saying  that 
Fermat's  principle  with  the  velocity  proposed  by  Rayleigh  is  wrong. 

He  see  however  that  the  rays  do  obey  Fermat's  principle,  the  correct 

expression  of  which  is  Eq.  36. 

IV.  THE  SAY  EQUATIONS  FOR  PLANAR  RAYS  IN  A  MEDIUM  WITH  v®  «  1. 

If  the  speed  of  the  medium  7  at  every  point  Is  much  smaller  than 
the  Local  speed  of  sound  c  (which  is  the  esse  for  most  propagation 
problems  in  either  the  ocean  or  the  atmosphere)  the  differential 
equations  for  the  rays  (Eqs,  19  or  28)  admit  considerable  simplifi¬ 
cations.  Keeping  only  first-order  terms  in  v  (as  compared  to  unity) 
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the  function*  N  and  ^ of  Eq,  19  reduce  to 

N  -  n;  X"  n(1  “ •  (39) 

We  ehell  conaider  a  further  simplification  by  assuming  that  n  and 
depend  only  on  x  and  y,  and  that  ^  Ilea  in  the  x-y  plane: 

n  -  n(x,y)  ;  ^ «  vE(*,y)ji  +  vy  <*.?)£  I  (40) 

where  ^  are  the  uaual  Cartesian  unit  vectors  in  the  x,y  directions. 
In  that  case  raya  which  start  out  in  the  x-y  plana  will  always  remain 
in  that  plane.  Such  rays,  therefore,  can  be  described  by  their  curva¬ 
ture  k,  given  by 

a'  ■  *£  •  (4i) 

where  the  unit  normal  vector, 

.  (42) 

which  is  perpendicular  to  the  unit  tangent  vector 


(43) 


hi  n«  uau.nvu  in  &a.  (ai  l  run  M  a  j  n .. 

F  "I-  >  UlSUUl.  xy 


from  Eq.  19a,  The  term  (r'*^')^  in  that  equation  ia  of  the  order 
nv  Jfc  ’  *nd  therefore,  compared  with  the  first  term  ,  is  negligible. 
The  third  term  can  be  worked  out  to  give 


W*  -  ay&oj}. 


(44) 
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with 

Vt  •  nv,(l  “&'•%)  »  1  •  xfy  •  C45) 

The  la at  tana  Q£n) ,  and  part  of  the  firat  term  ^r'dn/ds)  can  be  com¬ 
bined  by  uaing  the  idantitiea  dn/ds  •  and  x/S  +  y  /8  ■  1  to 

yield 

Jgn  -  (dn/daj^'  ■  (x'dn/8y  -  y'8n/dx)jS  .  (4b) 

Equation  19a  then  becomes 

ry^  -  (x'an/ay  -  y'WSx  +  8Vy/5x  -  avx/8y)£  ,  (47) 

whereupon  comparison  with  Bq.  41  gives  for  the  curvature 

k  -  n_1(x 'in/ay  -  y'an/Sx  +  3Vy/3x  -  3Vx/&y)  .  (48) 

For  ^  ■  0  this  equation  agrees  with  the  result  first  derived  in  Ref.  IS. 
Equation  48  is  still  quite  complicated.  It  is  a  second-order  differen¬ 
tial  equation  for  x(s),  y(s),  which  together  with  x'*  +  y,#  ■  1  can  be 
integrated  to  obtain  the  ray  paths  in  parametric  form  with  the  arc 
length  a  as  the  paraseter.  Or  one  can  eliminate  the  parameter  s,  and 
obtain  the  equation  of  the  ray  path  in  the  form  y  -  y(x)  directly, 
as  the  solution  of  a  single  second-order  equation  for  y.  In  either 
case,  however,  the  solutions  usually  would  have  to  bs  obtained  by  a 
numerical  integration* 

If  the  medium  is  simplified  somewhat  more  (as  .n  Ref.  1),  analytical 
solutions  of  Eq.  48  are  possible,  at  least  in  terms  of  a  quadrature. 
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If  the  index  of  retraction  n  and  the  velocity  of  the  medium^  are 
functions  of  only  one  coordinate,  say  y,  and  furthermore,  if 

«  v(y)^  has  only  a  component  in  the  x-direction,  then  Eq.  48  simpli¬ 
fies  further: 

K  *  n~1(x/dn/dy  -  dV/3y)  ,  (49) 

where 

V  -  nv(l  -  vx')  .  (50) 

Substituting  Eq.  50  into  Eq,  49,  and  using  the  relation 
(dn/dy)/n  -  -  (dc/dy)/c  one  can  obtain  the  following  expression  for 
the  curvature: 

h  -  c"l(dc/dy)(v-x')  +  (dv/dy)(2x'v~l)  ,  (51) 

which  is  identical17  to  Eq,  24  of  Ref,  1.  Thus,  even  though  the 
general  ray  differential  equation  in  Ref.  1  iB  incorrect,  the  small- 
velocity  approximation  is  correct.  In  particular,  the  analytic 
solution  in  terms  of  a  quadrature  fcr  a  medium  with  a  constant  wind 
(d(vc)/ay  »  0),  as  well  as  the  two  applications  (atmosphere  with  a 
constant  speed-of-aound  gradient  and  atmosphere  with  a  constant  tem¬ 
perature  gradient)  are  coirect.  The  reason  that  the  sections  of  Ref.  1 
which  deal  with  the  small-velocity  approximation  are  correct  is,  of 
course  the  fact,  that  the  term  in  Eq.  l»a,  whose  omission  is 

the  error  of  Ref,  1,  is  negligible  in  this  approximation. 
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APPENDIX  A 


In  this  appendix  a  brief  account  is  given  of  the  theory  of  the 
eikonal  equation  and  of  the  geometry  of  the  rays  and  wave  normals  for 
subsonic  as  well  as  supersonic  flows. 

The  wave  surfaces  are  described  by  tha  equation 

*  -  c0t  ,  ’  (Al) 

where  t  la  s  solution  of  the  eikonal  equation,  Eq.  I,  and  t  is  the  time 
from  an  arbitrary  origin.  From  this  it  easily  follows  that  the  unit 
vector  £  points  in  the  direction  in  which  the  wave  moves  and  that  p 
is  the  reciprocal  wave  speed  multiplied  by  c0 ,  Let  us  introduce  the 
following  notation  for  the  non-dimensional  wave  speed 

w  -  •  S  .  (A2) 

cp  ? 

It  is  a  quantity  which  by  definition  is  nonnegative.  The  eikonal 
equation  implies  that 

w  »  *  l  (A3) 

vhich  is  equivalent  to  Eq.  2.  This  equation  gives  the  speed  at  which 
a  wavelet  moves  provided  its  normal  and  direction  of  motion  are  known. 
One  could  use  it  in  principle  to  calculate  the  motion  of  a  wavefront 
if  it  were  not  for  the  ambiguity  in  the  sign  of  the  second  term  in  the 
right-hand  member  and  the  sense  of  the  unit  vector  o„ 
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To  clear  up  these  ambiguities  let  us  investigate  the  Cauchy  problem 
of  continuing  a  save  given  at  an  instant  of  time,  say  t  ■  0.  It  is 
convenient  to  introduce  a  local  right-handed  Cartesian  coordinate  sys¬ 
tem  at  the  wave.  Let  n  (not  to  be  confused  with  the  index  of  refraction) 
be  the  coordinate  normal  to  the  surface,  so  that 

v.  »  v»n  ^  0 

where  ^  is  the  unit  vector  normal  to  the  wave  and  pointing  into  the  half 
space  Into  which  points .  The  eikonal  equation  becomes 

5*  -  (l  -  i  v,  .  <A4) 


This  quadratic  equation  for  }s  has  the  solutions 


1 


1 


vB  ±  1 


(A5) 


Therefore,  provided  v,  •  1  ^  0  on  the  initial  surface,  there  are  two 
possible  solutions  to  the  Cauchy  problem.  One  may  Introduce  the  con¬ 
cept  of  "advancing"  and  "receding"  wavelets  with  respect  to  the  flow. 
The  plus  sign  refers  to  the  advancing  wavelet  and  the  minus  sign  to 
the  receding  wavelet.  Equations  AS,  A2  and  the  definition  of  p  imply 
that  the  wave  speed  is  given  by 


Kill 


(A6) 


This  implies  that  the  advancing  wavelet  has  a  greater  wave  speed  than 
the  receding  wavelet.  When  v„  “  0  there  is  no  distinction  between 
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advancing  and  receding  wavelets  and  the  local  motion  of  the  wavelet  is 
like  that  in  a  medium  at  rest.  Leaving  aside  for  the  moment  the  pos¬ 
sibility  that  v,  -  1  ■  0  on  the  initial  wave  surface,  one  may  now 
determine  the  sign  to  be  taken  in  Eq.  A3  or  equivalently  Eq.  2.  The 
results  up  to  this  point  are  given  in  Table  A1  as  the  first  four  entries 
for  each  of  the  four  possible  cases  (advancing  or  receding  wavelets  in 
subsonic  or  supersonic  flow).  It  la  clear  that  va  >  1  can  occur  only 
in  a  supersonic  flow.  When  gfjZ  >  0,  the  fluid  velocity  vector  points 
into  the  half  apace  into  which  the  wavelet  is  moving,  when  v^  <  0  the 
converse  holds. 

We  have  now  determined  what  the  sign  should  be  in  the  elkonal 
equation,  Eq.  2,  and  the  equivalent  Bq,  A3,  which  gives  the  wave 
speed  as  a  function  of  the  wave  normal  j^.  We  have  also  determined 
what  the  sense  of  the  wave  normal  p  is  In  relation  to  the  velocity  v. 

jm  *** 
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Instead  of  advancing  the  wave  normal  to  itself  using  Eq.  A3,  it  is 
more  convenient  to  advance  each  wavelet  along  the  ray  using  the  system 
of  ordinary  differential  equations  23a,  23b.  We  wish  to  study  these 
equations  in  greater  detail  in  order  to  remove  the  ambiguity  in  the 
sign  of  Eqs.  24a  and  24b  (see  Refs.  11  and  8),  and  also  to  find  a 
simple  geometric  construction  relating  rays,  wave  normals,  ray  speeds 
and  wave  speeds. 

To  the  characteristic  equations  23a,  23b  one  may,  using  Eq,  21c, 
add  the  following  relation: 

l 

M  •  2n(n  -  v,pt)  •  ±2n(p,p,)*  -  ±2op  ,  (A7) 

a  t 

where  use  has  bean  made  of  the  elkonal  equation.  This  equation  com" 
bined  with  Eq.  Al  gives 

co  7^  "  “  *2r>P  ,  (A8) 


or 


dt 

c  dT 


(AS) 


which,  apart  from  the  ambiguity  in  the  sign,  relates  the  parameter  t 
to  physically  significant  quantities.  The  rule  of  signs  to  be  fol~ 
lowed  is  the  same  as  in  the  eikonal  equation  written  in  the  form  of 
Eq.  A3. 

One  may  now  east  the  characteristic  equations  23a  in  a  geometri¬ 
cally  perspicuous  form.  Introducing  time  as  independent  variable 
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(A10) 


instead  of  t  one  gets  from  Eq.  23a 


i_  fill.  -  i-  ii 

ce  dt  Cq  dt 


■Mr*  ’•} 


or 


I 

c  dt 


■  v,  ±  • 

P 


(All) 


In  the  notation  of  this  paper  this  may  be  written 

ss/  ■  & ±  jL  •  (A12) 

This  Is  the  generalization  of  Eq.  3  and  applies  also  to  supersonic 
flow.  The  results  obtained  are  collected  in  Table  Al.  This  permits 
us  to  give  the  relationship  between  ray,  wave  normal,  ray  speed  and 
wave  speed  in  graphical  form  (see  Figures  Al  and  A2  for  the  subsonic 
and  supersonic  cases  respectively).  In  either  case  the  circle  is  of 
unit  radius  centered  at  0  and  the  line  A$  represents  the  velocity 
vector 

In  the  subsoni'.  case,  (v  <0),  which  is  normal  in  acoustics,  the 

•* 

situation  is  quite  simple  as  indicated  in  Fig.  Al.  The  vector  A3 


represents  the  ray.  Its  langth  Is  q,  the  dimensionless  ray  speed. 
The  ray  is  made  up  of  the  vector  sum  of  the  velocity  £  and  the  unit 
wave  normal  ^represented  by  OB  (see  Table  Al).  The  projection  of 
%  on  plus  unity  is  the  wave  speed  represented  by  the  length  of  CB. 


(Mt  a  w/f  kh  Ual  f  aLa  ^  i  £ _ _ _ a_l J s  _  -  _  * 

*»*v  a  *9ttw  iimaa  w a.  baas  WAAiiUUtiKfc  W1VC  UA  LUC  UllCie  AepEeBORlB  ftQV&nC' 

ing  wavelets  and  the  left  half  receding  wavelets. 

The  relationships  for  the  supersonic  csss  (v  >  1)  are  more 
involved  and  aru  indicated  in  Fig,  A2.  A  line  drawn  from  A  to  a  point 


23 


on  the  circle  represents  a  ray  as  before.  One  has,  for  example,  AB 
and  aB '  two  rays  with  the  same  direction,  one  fast  and  one  slow. 

The  unit  normal  £  pertaining  to  A?  is  61?  and  that  pertaining  to  AB ' 
is  The  wave  velocity  vector  wg^  pertaining  to  the  two  rays  are 

CB  and  c'a  respectively.  Ve  see  that  by  considering  only  the  plus 
sign  in  Eq.  3  we  confine  ourselves  to  the  fast  rays  in  the  supersonic 
case.  In  the  supersonic  case  we  note  further  that  the  angle  3  between 
the  ray  and  the  velocity  vector  must  satisfy 


where 


lei  *  u 


sin  4  ■  - 
v 


(A13) 


and  i*  is  called  the  Mach  angle.  Thus  one  may  not  start  a  ray  with 
an  arbitrary  slope. 

From  the  graphic  representation  one  may  derive  various  relation¬ 
ships  by  trigonometry.  Ve  give  sane  below. 

The  non-dimensionel  ray  speed  as  a  function  of  the  angle  the 
ray  makes  with  the  velocity  vector  is 


A 

q  ■  v  cosB  ±  v  -  oin*!y  . 


(A14) 


When  v  <  1  (subsonic  flow)  the  plus  sign  is  to  be  used.  When  v  >  1 
(supersonic  flow)  the  two  signs  give  the  fast  and  alow  rays  respec¬ 
tively. 
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The  cosine  of  the  angle  between  ray  and  wave  normal  (it  is  always 


an  acute  angle)  Is  given  by 


which  is  Identical  to  the  quantity  S  defined  before, 
less  wave  speed  may  be  found  from 

w  -  <£,'.£  . 


(A15) 

The  dimension- 

(A16) 


In  conclusion  one  may  now  see  what  the  meaning  of  v,  •  1  •  0  is, 
a  case  that  was  excluded  from  the  discussion.  It  pertains  to  a 
receding  wavelet  carried  along  by  the  fluid  so  that  its  normal  speed 

-f 

vanishes.  It  corresponds  to  the  rsy  AB  in  F%.  AZ.  Clearly  on  such 
a  wavelet  (grad  f j  •  fl  ■  •  and  the  ray  AB*  is  tangent  to  the  wavelet. 
At  such  a  wavelet  Eq.  A1  cannot  hold.  In  mathematical  language  the 
•ikonal  function  I  must  be  singular  there.  The  surfaces  on  which 
v  -  0  appear  as  branch  surfaces  of  4,  joining  the  fast  and  alow  waves. 
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TABLE  Ala  The  basic  equations  for  advancing  and  receding 
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fig,  A-2:  Relationship  between  lays  and  wave  normals  at  a 
supersonic  point,  A.  W, :  advancing  wavelet, 

R.  W. :  receding  wavelet. 
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Two  different  foxiuulations  axe  presented  in  terns  of  ordinary  second-order 
vector  differential  equations  for  acoustic  rays  in  a  three-dimensional  medium 
moving  with  an  arbitrary  velocity  (either  subsonic  or  supersonic),  and  having 
an  arbitrary  index  of  refraction.  The  first  formulation  has  the  arc  length 
of  the  ray  as  independent  variable,  whila  the  second  one  is  given  in  temp 
of  a  canonical  varible  and  is  equivalent  to  Keller's  Hamiltonian  formulation, 

(J.  B.  Keller,  J.  Appl.  Ffcys.  2£.  938-5*7  (1954)).  >v 
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